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Abstract 

We consider the propagation of massive-particle de Broglie waves in 
a static, isotropic metric in general relativity. We demonstrate the exis- 
tence of an index of refraction that governs the waves and that has all 
the properties of a classical index of refraction. We confirm our inter- 
pretation with a WKB solution of the general-relativistic Klein-Gordon 
equation. Finally, we make some observations on the significance of the 
optical action. 



I. INTRODUCTION 

In the case of static, isotropic metrics, the optical-mechanical analogy leads to 
a remarkable simplification of the geodesic equation of motion. The trajectories 
of both massive and massless particles can be parameterized by one and the 
same parameter, the optical action A. Moreover, for both photons and massive 
particles, the equation of motion assumes the form of Newton's second law. Thus 
the exact general-relativistic problem becomes formally identical to the classical- 
mechanical or, if one prefers, to the classical geometrical-optical, problem. 

This formalism also provides an easy way to introduce quantum-mechanical 
wave packets into the context of general relativity. In this paper we consider 
the propagation of massive-particle de Broglie waves in a gravitational field 
represented by a static, isotropic metric. We shall use simple, but rigorous, semi- 
classical arguments to demonstrate the existence of an index of refraction that 
governs the waves. Our expression for the index of refraction will be confirmed 
by a WKB solution of the general-relativistic Klein-Gordon equation. We shall 
see that the same index of refraction may be used in both the classical-particle 
and wave-mechanical pictures of propagation. Finally, this treatment will shed 
new light on the physical significance of the optical action. 

A. The gravitational field as an optical medium 

Our point of departure will be an exact formulation of the geodesic equation of 
motion which is nevertheless of Newtonian form. The remainder of the Intro- 
duction will summarize this formulation of the geodesic problem. This summary 
will serve not only to fix the notation but also to introduce readers to the opti- 
cal action, the physical quantity that plays the central role in all that follows. 
Throughout this paper, we suppose that the metric is static and isotropic, so 
that the line element is of the form 

ds^ ^ n\r)cldt^ - '^-\r)\dr\\ (1) 

where ft and <& are functions of the spatial coordinates r = (r, 0, (j)) or (x, y, z) 

and where cq is the vacuum speed of light. Many metrics of physical interest 

can be put into this form, including the Schwarzschild metric. 0] The isotropic 

coordinate speed of light c(r) is determined by the condition that the geodesies 

be null {ds = 0) : 

dv 
c(r)- - == co$(r)17(r). (2) 

Thus the effective index of refraction for light in the gravitational field is[0 

n(r) = ^-^n-^. (3) 



B. Relativistic geometrical optics in Newtonian form 

This index of refraction may be used in any formulation of geometrical optics. 
For example, the shapes of light rays are governed by Fermat's principle: 

S f \dl = 0. (4) 

6 represents a variation of the integral produced by varying the path between 
two fixed points in space, Xi and X2. dl = \dr\ is the element of the path of 
integration in the three-dimensional space. From Eq. (4) one may derive other 
formulations of geometrical optics, including the "-F = TTia" formulation. In 
this formulation, the equation of the ray assumes the form of Newton's law of 
motion: 

^=v(l„"ci). (l.ght) (5) 

This is an exact formulation of geometrical optics that was originally applied 

to classical isotropic media, [pi However it is equally valid for static, isotropic 

metrics in general relativity.jH] Eq. (5) reduces relativistic geometrical optics 

to the form of classical Newtonian mechanics. The left side of Eq. (5) has the 

form of an acceleration. The right side has the form of a force and the effective 

potential-energy function is — n^CQ/2. The key feature of this formulation of 

geometrical optics is the use of the optical action A as independent variable. A 

is defined by 

dr 

— =ncQ. (light) (6) 

dA is related in a simple way to the coordinate time dt elapsed as light moves 
an infinitesimal distance dl along the ray: since dA — dl/ncQ and dl ~ codt/n, 

dA = dt/n^. (7) 

C. Relativistic particle mechanics in Newtonian form 

The orbits of massive particles are obtained by requiring that they be geodesies: 

(•X2,t2 

S ds = 0. (8) 

Jxi.ti 

Here S represents a variation of the integral produced by varying the path be- 
tween two fixed points in spacetime, (xi,ii) and (x2,i2)- We wish to cast 
this variational principle into the form of Eq. (4) in order to make the optical- 
mechanical analogy as transparent as possible.]^, || 0] With the use of Eq. (1) 
and (3), Eq. (8) can be written in a form analogous to Hamilton's principle: 



■I 



t2 

L{x„Vi)dt^O, (9) 



where the effective Lagrangian is 

L{x,,Vi) = -mcl^ [1 - V^n'/cl] ^'^ . (10) 

In the Lagrangian, 51 and n arc functions of the spatial coordinates alone, V^ = 
X]i=i ^i if ^^ choose to work in Cartesian coordinates, and Vi denotes dxi/dt. 
The factor of the rest mass m has been included for dimensional convenience. 
The canonical momenta are 

p,=mnn^[l-V^n'^/clY^''^V^. (11) 

The effective Hamiltonian is 

H = mcln[l-V^n^/clY^'\ (12) 

or, if expressed in terms of the momenta, 

H = mcl [n" + p'/n'm'c^^ ^'^ . (13) 

77 is a constant of the motion, whose value may be calculated from Eq. (12) 
if the particle coordinate speed V is known at one point on the path. A little 
algebra applied to Eq. (12) gives the speed function V{y): 

y(r) = con-i[l-m2c4r!2/H']'^'. (14) 

Eq. (14), for massive particles, corresponds to the familiar v = con~^ for light. 
From Hamilton's principle, Eq. (9), one may derive Jacobi's form of Mau- 
pertuis's principle: 

/•X2 

5J pdl^O, (15) 

where now the path of integration though three-dimensional space is varied, 
subject to conservation of energy, between two fixed points in space, xi and 
X2, but the times at the end points need not be held fixed. From Eq. (11) and 
Eq. (12) we have 

p, - Hn^Vjcl (16) 

and thus Eq. (15) can be written 

sT^dl^O. (17) 

Here n^V is to be regarded as a function of the spatial coordinates. 

Eq. (17), which governs the orbits of massive particles in static, isotropic 
metrics, is formally analogous to Maupertuis's principle, which is a suitable 
foundation for the Newtonian mechanics of a point-particle in a static, isotropic 
potential. But Eq. (17) is also analogous to Fermat's principle, Eq. (4), which 
governs the shape of light rays in isotropic media. Thus, corresponding to 



any formulation of geometrical optics that is derivable from Eq. (4), there is a 
formulation of relativistic particle mechanics derivable from Eq. (17). 

Indeed, the exact general-relativistic equation of motion for a massive par- 
ticle may be written 

(fr /I \ 

-— ^ = V ( -n'^V'^ I . (massive particle) (18) 

This equation is precisely of the form of Newton's second law: acceleration = 
negative gradient of potential energy. Note that, because V = cq/ti for light, 
Eq. (18) reduces to Eq. (5) as a special case. 

The independent variable A playing the role of the time in Eq. (18) is the 
optical action, defined by 

dr 

v?V. (massive particle) (19) 



dA 

But since V = |dr/rfi| it follows that, along the trajectory, 

dA = dtjn^. (20) 

Thus the same independent variable A is used to parameterize the path, whether 
the particle be massive or massless. The only difference in the treatment of 
massive and massless particles is in the form of the potential energy function, 
C/(r) = — n^Co/2 or — n^y^/2, which appears under the gradient operator in 
Eq. (5) or Eq. (18). 

II. DE BROGLIE WAVES IN THE 
GRAVITATIONAL MEDIUM 

With these results in hand, we are ready to introduce the de Broglie waves of a 
massive particle into the medium that represents a static, gravitational field. Let 
these waves have phase speed Wp(r), wave number fc(r) and angular frequency w, 
related by v^ = tj/fc. Because of the similarity of Eq. (4) and Eq. (17) (as well 
as the similarity of Eq. (5) and Eq.(18), we might speculate that TV = n^F/co 
plays the role of an index of refraction for the massive particles. We shall see 
that this is indeed the case. 

In the geometrical-optics limit, it is possible to speak of a ray, which is the 
trajectory of the point-particle. The ray will be the path satisfying Fermat's 
principle, which we write in the form: 

(5 / kdl = 0. (21) 

And now we may reason in exactly the same way as Louis de Broglie in his thesis 
of 1924.11 The paths predicted by Eq. (17) and Eq. (21) will be the same if the 
two integrands are the same functions of the spatial coordinates. However, they 
may differ by a multiplicative factor /. Thus we require 

Hn^V/cl = fk. (22) 



A. The de Broglie wavelength of a massive particle 

To determine / let us pass to the empty-space lim.it in which H — > ?7iCg7 = 
mcl{l - V'^/cl)-'^/'^ and n^ -^ 1. Eq. (22) becomes mV-/ = fk. This will be 
precisely the de Broglie relation of a relativistic free particle ii f = h = h/2TT, 
where h is Planck's constant. Thus we have 

Hri^V/cl = fik, (23) 

which is just what we would have obtained if we had simply put 

p^nk (24) 

using Eq. (16). 

Since {H/cq^ p) and [lo / cq, k) must each transform as a four-vector, we must 
also have 

H = huo. (25) 

Thus, given a wave of angular frequency a;, the frequency remains constant 
along the whole trajectory or ray. For a particle in an unbound orbit, to is 
everywhere equal to the value of w at spatial infinity. The optics of de Broglie 
waves in a static, isotropic metric is therefore analogous to the wave optics of 
classical isotropic media. Of course, it must be emphasized that w is a coordinate 
frequency^ and not a directly measurable frequency, m is 27r times the number 
of wave crests crossing through a fixed point r per unit coordinate time t. From 
Eq. (|]), the frequency that would be measured by an observer at r is a) = uj/Vtir). 
From Eq. (P3|), the wavelength of the massive particle in the gravitational 
field is 

A=^- (26) 

It should be kept in mind that this is the coordinate length of the wave only. The 
physically measurable wavelength is obtained by applying the metric of Eq. (1). 
Thus the measurable wavelength is A = A/<i>. 

Using Eq. (14), we can cast Eq. (24) into the form 

^- m' (27) 

which specifies A entirely in terms of the space functions n and VL and the 
constant of the motion H . 

B. Index of refraction for massive-particle de Broglie waves 

From Eq. (p7|) we see that the quantity 



A^l-^-^ (28) 



is a constant everywhere in the medium. Now, Eq. (Eq), which governs the 
de Broglie waves of the massive particle, is analogous to the familiar rule from 
ordinary light optics (An — constant). Thus the wave-optics rule for massive 
particles may be expressed 

A7V = constant, (29) 

where the index of refraction iV for the de Broglie waves of massive particles in 
the gravitational field is given by 



„ = „yi_!2^. (30) 

Eq. (pQ) is an exact result, applying to fields of arbitrary strength. The gravi- 
tational wave-optics of massive and massless particles differ only in a dispersion 
effect, since the index of refraction for massive particles depends on the param- 
eter m/H. From Eq. (f2), note that for hght ni/H -^ and then N ^ n. With 
the use of Eq. (14), Eq. (|3^) may also be written 

N = nV/co, (31) 

as we expected. With use of Eq. (16), we also have 

N = cap/H. (32) 



Finally, with use of Eq. (25), the index of refraction may be written explicitly 



as a function of position and frequency: 



N{r,u;)=n{r)Jl ^f-^i. (33) 

By virtue of Eq. (p3l|), the geodesic equation of motion of the massive particle, 
Eq. (18), may be written 

(fr (\ \ 

—— = V I -A^^Cg I , (light and particles) (34) 

and the optical action A may be defined by 
dv 



dA 



Nco- (light and particles) (35) 



These relations are equivalent to Eq. (18) and Eq. (19); but now, remarkably, 
they take on the forms of Eq. (5) and (6). All these equations are, it should be 
stressed, exact general-relativistic formulations in the class of metric under con- 
sideration. In the special case of static, isotropic potentials or media, Eq. (p4|) 



and (35) therefore cover all of the following: classical geometrical optics, clas- 
sical Newtonian mechanics, general relativistic geometrical optics and general 
relativistic particle mechanics. The effective potential energy function is always 
just the negative square of the index of refraction. 



C. Phase and group velocity 



The phase velocity of the waves is Vp — uj/k = H/p, or, with use of Eq. (32), 

vp = co/N, (36) 

a relation which demonstrates once again that N is playing its proper role as 
an index of refraction. Inserting p ~ Tik and H — hui into the square of Eq. (n3f) 
we obtain 

ri^u^ = m,^4n^ + cln^k^/n^. (37) 

Differentiating with respect to k, we obtain 



■-0 



k dk -n?' 



(38) 



VpVg = cl/ri^, (39) 

where Vg = du/dk is the group velocity of the de Broglie waves. Using Eq. (p6|), 
we find 

Vg = c^N/n^. (40) 

Substituting the expression for N from Eq. (pi]), we see that 

Vg = V, (41) 

the coordinate velocity of the classical point-particle. Using Eq. (20), we can 
also write Eq. (5^) as 



VpVg — Cg 



dA 
dt' 



(42) 



D. Red shifts of massive particles 



Another way to see that N really is playing the role of an index of refraction 
for massive particles is to show that it yields the correct gravitational red shifts. 
Let a particle of mass m travel in free fall between position ri and position r2. 
The coordinate wavelength A associated with the particle varies according to 

iV(ri)A(ri) = iV(r2)A(r2). (43) 

The physical or metric length of the waves is A = A/$, so the rule governing 
the redshifts of massive particles in static metrics is 

iV(ri)$(ri)A(ri) = A^(r2)$(ri)A(r2), (44) 

or 



A(ri)^r!-2(ri) - 4/H^ = A(r2)^r!-2(r2) - 4/H^. (45) 

Eq. (45) and (46) are entirely analogous to the rules governing the redshifts of 
light (for which n stands in place of iV).p| 



III. WKB SOLUTION TO THE GENERAL 
RELATIVISTIC KLEIN-GORDON EQUATION 

The index of refraction for massive particles given by Eq. (pO|), (|3l| ) or ( p2|) 
has been obtained from semi-classical considerations. We would like to show 
that the same expression results from examination of the appropriate quantum- 
mechanical wave equation. The Klein-Gordon equation for a spinless massive- 
particle wavefunction in flat space results from making the substitution p^ -^ 
—ihd^ in the Minkowski-space mass-shell constraint g^'^p^iPu — (mco)'^ — and 
allowing the resulting operators to act on a scalar wavefunction ^. In the 
general-relativistic generalization of this wave equation to curved spacetime, 
the momentum operator must be replaced by the covariant derivative, so that 
Pf^ —f —iKS/f^. Making this generalization, we arrive at the general relativistic 
Klein-Gordan equation (GRKGE), a good candidate for a wave equation for 
massive particles of zero spin:pO[ 

(^.g'^'^V^V, + ^^V(a;) = 0, (46) 

where il'i^) is our scalar field, 5^" is an element of the metric tensor, and we 
follow the usual convention of implied summation over repeated indices. 

Now, the covariant derivative operating on a scalar function just gives the 
usual partial derivatives: VyV' = dv'ip- However, di^tp forms a covariant vec- 
tor. And when a covariant derivative operator acts on the covariant component 
Qu{x) of a vector, we have 

V^Q,(.t) = d^Q,{x) - r^,(x)QA(x), (47) 

where the F^^ are the connection coefficients defined by 

r^. = ^9^" {d^g.. + 9.5p^ - 9«5m.) • (48) 

Expanding out the double covariant derivatives, we may write the GRKGE 
in the more explicit form 

[g^''d^d,~g^^''Tl,dx + {mcoinf]i;{x) = 0. (49) 

This equation is the jumping-off place for most treatments of scalar matter 
waves in curved spacetimes.pT[ 

In the static, isotropic metric of Eq. (1), g^y = diag.(ri^Co, — $^^, — <i>^^, -$^^) 
Thus the inverse tensor g^^" = diag.(r2~^C(^^, -<I>^, — 4>^, — $^). And g = det g^^ = 

Thus, with the aid of the identity |12] 



'T;:. = - i-9)-"'d, li-gf'^g^"] , (50) 



we obtain 

where we define 
and 



n2 av „2 



df^ 



V^V - VC • VV' + fee (r) V' = 0, (51) 

f (r) = ln(f^/$) (52) 



fcc (r) = mco/h^. (53) 

h/mcQ is the invariant, or physical, Compton wavelength of the particle. From 
the considerations in Sec. II. D., $/i/mco may be interpreted as a coordinate 
Compton wavelength. Thus kc (r) is the coordinate Compton wave number. 



We can remove the term V^ ■ V'i/' from Eq. (51) by putting 

i;{r,t) = f{r)<j>{r,t), (54) 

where 

/(r) - i^/nf' = e-^l\ (55) 

Then (r, i) satisfies the differential equation 

cl 9t2 
where we define 



V2(/)+ [/c2(r) + ,^(r)]0 = O, (56) 



,y(r)^iv2e(r) + i|Ve(r)|2. (57) 

Note that Eq. ( pq ) is essentially the flat-spacetime wave equation for a scalar 
particle in a spatially varying index of refraction n and with a spatially varying 
Compton wave number k^. However, there is an additional term 77(r) in the wave 
equation which results from the geometric cross term g^^^\y in the product of 
the covariant derivatives. In the weak field limit this term is typically dropped 
[O when seeking an approximate solution to the wave equation. However, it is 
not much extra trouble to carry this term along. In the next section we develop 
a WKB solution to Eq. (Im 



A. The WKB expansion 

Let us look for a stationary solution to Eq. ( pq ) of the form 

i 



>{r,t) — Cexp 



/(r,i) 



(58) 



with C constant and where we will take the phase to be of the form 

5(r,t) = 5(r)-iJi, (59) 



where H is constant. Substitution of Eq. ( |5§| ) and Eq. (|59|) into Eq. ( |5q ) shows 
that we will indeed have a solution, provided that 

p^ -{\7Sf+in\7^S-h'^r]{r) = 0, (60) 

10 



where we have used Eq. (13) and Eq. (3). 

We now expand S'(r) in a power series in h: 



5(r) = ^?."54r). (61) 



n=0 



Substituting this into Eq. (60) yields the set of equations 

0(?i") : {VSof-p^^ ^ (62) 

0{h^) : 2VSo ■ VS"! = iV^So (63) 

0{h'^): 2VSo-WS2 = iV^Si-{VSif-T]{r) (64) 

n-l 

0(?i"2^3) : 2VS'o-VS'„ = iV2S'„_i-^VS'j- VS'„_j. (65) 

B. The eikonal equation 

Eq. (|62|), the eikonal equation, governs the dominant contribution to the phase 
of the wave. By direct integration we have 

So{v)^ J p(r).dx, (66) 

where 




p(r)^|p| = -n^/l-(^^) . (67) 

This form for p comes from rearrangement of Eq. (13). 

The surface Sq (r) = constant is a wavefront with normal in the direction of 

VSoir) = p(r). (68) 

Following Holmes JMJ, let us characterize a wavefront by coordinates r = 
5c{l, a, (3), where I is the arc length along the ray, normal to surfaces of constant 
5*0, and a and /3 are coordinates used to parameterize the wavefront surfaces 
(e.g., spherical coordinates). dx./dl is the unit vector that is tangent to the ray 
or trajectory. From Eq. (16) we see that dx/dl is in the direction of p(r) and 
hence, by Eq. (|68|), in the direction of VS'o(r). Thus 

dx ^ pW ^ V5o(r) 

dl p |V5o(r)|- ^ ^ 

For any F{r) we can write 

dF dx „^ .„„. 

Putting Sq for F in Eq. (JTy) and using Eq. ( p^ we have 

^^f.V5..|V5„^ (71) 

11 



And thus, by Eq. (|68|), 



^o(/,«,/3) 



pdl, 



(72) 



which is another way to write Eq. (|6^) • J pdl is proportional to the optical path 
length along a segment dl of the ray. And thus p is proportional to the index of 
refraction. |15|| So here we have a more detailed justification of Eq. (pa). 



C. The transport equation 

We now need to solve the transport equation, Eq. (|63|). Putting Si for F in 
Eq. (|ol), we have 



dl 
where we have used Eq. (p&l). Then, by substitution, Eq. ( |63| ) becomes 

dSi _ i ^2c. 



with solution 



y'So 



(73) 

(74) 
(75) 

(76) 

d{l,a,(3) ^^^^ 

is the Jacobian of the transformation from the curvilinear ray coordinatesf/, a, /3) 
to Cartesian coordinates. Then, upon substitution of Eq. (u&) into Eq. dTq), one 
can perform the integration to obtain 



Si{l,a,(3)='- f dl — 
2 7 P 

A short calculation shown in the Appendix reveals that 

W^So d{pj)/dl 



where 



J = 



pj 
9x 



^' <*"■«' = 5'° (r^)4''<"' 



(78) 



where ^(r) = In ( pfp )j(r ) ) ^^^ '"o = ^\i=o- To the lowest order (0(?i°) in the 



{ro)J(ro) 

phase and amplitude), we have found the WKB approximate solution 



$(r) p(ro)J(ro) 



^n p(r).rfx-i/t 



(79) 



where the subscript 1 on tpi (r, t) indicates that we have carried out the WKB 
expansion to Sn=i- 

Note that Eq. d79) takes into account, to lowest order, the term ^(r) = 
ln(ri/$), which arises from the covariant derivative term .9^T^^. What has 



12 



been left out are higher order terms involving V^^ (r) and V^ (r) in 77 (r), which 
come from the quantum corrections terms, Sn>2 (which we deal with next). Note 
also that this solution is valid for strong gravitational fields. [O But truncating 
the solution at 5*1 does impose the restriction that p^ 3> h^rj (r), or 



A|VC|, A|yv^|«l. (80) 

Since f (r) is slowly varying, the above conditions can hold even reasonably close 
to the horizon of a black hole. We investigate this in the next section. 

D. Estimation of terms in Schwarzschild metric 

If we had defined (j){r,t) = u (r) exp{—iHt/h) our wave equation, Eq. (p6|), 
would have the form of the Helmholtz equation: 



V2w(r) + fc2(r) 



^ vi^) 



u(r) = (81) 



where p{r) = h fc(r), fc(r) = fco A^, with fco = i^o/cq = H/Ticq. Therefore we want 
to consider the order of magnitude of the terms | V^j^/fc'^ and | V^^|/fc^. 
For the Schwarzschild metric in isotropic coordinates we have: 

n = 
$ = 



1-1/p 
1 + 1/p 

1 


(1 + 1/pf 
(1 + ypf 


(1 - i/p) 

ln(l - l/p') 



k = kony/1 - n^/H'^, (82) 

where iJ'^ ee H/{mcl) and p = r/r^, with r^ = GM/cl = 0.74 M/Mq km (half 
the Schwarzschild radius). 

For massless particles, H' — > 00, so that k — k^n and we obtain 



|V^(r)P _ 4 1 1 

fc'(r) {k^r^f P^ {1 + 1/ pf 

|V2e(r)| _ 2 1 (l + l/p2) 



fc2(r) {karsfp^{\ + llp) 



(83) 



These ratios are finite for all values 1 < p < 00 (the valid range of the isotropic 
scaled radius p) , even though the numerators and denominators of the left hand 
sides of Eq. ( |8^ ) each separately diverge as p -^ 1. From Eq. ([8^), k{p) -^ fco as 
p —>■ 00, so fco is the usual wave number of the particle at spatial infinity. For 
ordinary wavelengths, fcor^ ^ 1, since Vs ~ km, so the above ratios remain very 
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small even down to p = 1. Even if we were to consider ultra-long wavelengths 
such that koTg ~ 1, the ratios in Eq. (p3) could still be made small for values of 
/9 ~ 2, i.e. r = 2rs, the Schwarzschild radius. In this case, the term ^^(r) arising 
from the covariant derivative term f/'^'^F^^ is negligible compared to fc^(r) all 
the way down to the Schwarzschild radius, and for all intents and purposes, our 
wave equation for massless particles is of the form 



V^M(r) + /c^(r)u(r) = 0. 



(84) 



However, in the next section it causes no great difhculty to carry terms involving 
r](r) along formally. 



E. Quantum corrections 

The first occurrence of rj (r) appears in Eq. (pi). Using Eq. ( [70| ) we can write a 
first order equation for 6,82 /dl whose solution is given by 



S2{l,a,f3) = -J ^ 



ivV(r)-i(VMr))V^(r) 



(85) 



Since ^2 is purely real, this is an OJh ) correction to the phase. 

The remaining equations, Eq. (pa), can be formally solved to give 



S„{l,a,(3) ^ Y 






(86) 



where the terms in the brackets involve only the previously determined functions 
Si, S2, ■ ■ ■ , Sn-i- The first correction to the amplitude is 0{h ) and is given by 



/■' dl 
S3il,a,P) = ij — [v252-V/i(r)-V52] = iS'^,. 

Putting this all together, we have to 0{h ) in the phase and amplitude 



(87) 



V'3(r,i) 



C^ 



/£(r) / p(ro)J(ro) 
n{r)V p(r)J(r) 



e^p{~h^S',{r)) 



exp 



/ p(r) ■dx-Ht + h'^S2{r) 



(88) 



where the subscript 3 on tp^ (r, t) indicates that we have carried out the WKB 
expansion to Sn=3- 

The lowest-order contribution to the phase of the wavefunction in Eq. ( p8|) 
is just —h~ times the classical action 



(89) 



Sci = / Ptj.dx'^, 
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where p^ = {H/co,p). If this four- vector is substituted into the mass shell 
constraint g^'^p^Pu — (mco)^ — 0, we obtain the square of the Hamiltonian, 
Eq. (p^. Using Eq. ( p^ ) and the definition of T^ = dl/dt, we can also write 
the classical action as 5c/ (r) — H/cq /"" {N — cq/V) dl along the geodesic. For 
massless particles, N ^ n and V -^ co/n so that Sd — > 0. Thus, the quantum 
correction term 52 (r) in Eq. ( |88| ) become the lowest order contribution to the 
phase of the wavefunction in this case. 

In 1979, Stodolsky ||lg argued that the wave function of a spinless quantum 
mechanical particle in a arbitrary metric ds^ — g^,ydx^^dx'^ should take the 
form ipiv^t) = A exp[— z5ci(r,i)/?i], with Sd given by Eq. (|8S|). This expression 
for the phase of a quantum mechanical wavefunction in curved spacetime have 
been used recently to compute the effects of mass oscillations from supernova 
neutrinos |]1^ . Eq. (|8q ) reveals that the phase of the wavefunction for a massive 
spin particle only has the form given by Eq. (|8|) to lowest order in Ti. The 
Klein-Gordon equation, which is second order in space, cannot have a solution 
in which the phase is given exactly by Eq. (B9^, since two spatial derivatives of 
such a phase would produce terms like V • p, which are not present in the wave 
equation. It is interesting to note that the wave function for a spin 1/2 particle 
in a static metric has precisely the form proposed by Stodolsky, Eq. (p9[) |18{ . 
This result follows from the first-order nature of the Dirac equation in curved 
spacetime. 

IV. THE OPTICAL ACTION 

A. The optical action and the proper time 



Eq. (35) defines the optical action for both massive and massless particles. 
Eq. (7), which also is valid for both massive and massless particles, provides 
the relation between the optical action and the coordinate time elapsed along 
a trajectory. It will be convenient to state explicitly the connection between A 
and the most commonly used parameter for massive particle trajectories, the 
proper time r. From Eq. (|l]) and Eq. (12), the relation between the elapsed 
proper time dr(= ds/co) and the optical action dA may be written 

dA = ^ dr. (90) 

r is not a suitable parameter for light rays since, for massless particles, dr -^ 0. 
However, in this case H/m ^ cx) in such a way that the product {H/m)dT 
remains finite, as we see from Eq. (7). This is why the o ptical action may be 
used for parameterizing the trajectories of both massive and massless particles. 

B. The optical action in the WKB calculation 

The optical action dA is related very simply to the measure dl/p, which appears 
in the integrals for the higher order WKB phase terms in Eq. (fZSh, Eq. (p5[) 
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and Eq. (|8|). From the definition of A [Eq. (g|)], we have dl = NcodA. Using 

Eq. (^, we obtain 

— M 

dA^^^dl. (91) 

%P 

Thus, the phase in the expressions for 5„>i is integrated with respect to the 
optical action. 

C. The optical action as Born and Wolf's r' 

In their discussion of geometrical optics, Born and Wolf [|l9[ introduce an oper- 
ator 

T^ = V^o • V, (92) 

or 

where r' is a certain parameter that specifies position along a light ray. (Born 
and Wolf actually write r rather than r'. We add the ' to avoid confusing this 
parameter with the proper time in general relativity.) Applying this operator 
to the eikonal 5*0 itself gives 

f^-(V^o)^=P^ (93) 

and thus 

dSo=p^dT'. (94) 

However, from Eq. (|7^) we have 

dSo=pdl. (95) 

Comparing these two expressions, we see that 

dr' = -. (96) 

P 

From this it follows that the operator defined by Eq. ( p2^ is pdi, that is, p times 
the derivative along the ray. Comparing Eq. (^^ with Eq. (pi|), we see that, 
apart from a factor of H/cq, Born and Wolf's dr' is the optical action dA. 

D. The optical action as a wave-and-particle action 



Eq. ( |42| ) also provides some insight into the significance of the optical action. 
Multiplying Eq. ( ^ ) by dt we have 

dA = Vp {vgdt)/cl = Vp dl/cl, (97) 

or in words (ignoring the factor of c§), 

change in ^ = phase speed of waves x displacement of particle. 



16 



Alternatively, this may be written 

dA = Vg {vpdt)/cl, (98) 

change in ^ ~ speed of particle x displacement of phase wave. 

These expressions for dA bear a rather intriguing relation to the classical action 
of a Newtonian particle (= speed of particle x displacement of particle). The 
optical action A for a massive particle has, as it were, a foot in the wave regime 
and a foot in the particle regime. 



V. SUMMARY 

In a static, isotropic metric, it is possible to define an index of refraction not only 
for massless, but also for massive particles. The index of refraction N is given by 
Eq. (|3^) , ( pi] ) or (^2|) . This index of refraction works as it should for de Broglie 
waves in the physical-optics rules of Eq. (|2^) and Eq. ( JSq ) . Thus one may work 
out physical-optics applications, including interference and diffraction, using 
this A^.p6| Passing over from physical optics to geometrical optics, we attain 
the point-particle regime. In the point-particle regime, the trajectory of the 
particle (whether massive or massless) is expressed by an equation of motion 
formally identical to Newton's second law, Eq. (p3), in which the independent 
variable is the optical action, defined by Eq. (pa)7 and in which the potential 
energy function is — A^^Cq/2. We have confirmed this index of refraction, and 
explored the limits of its validity, by giving an explicit WKB solution to the 
general-relativistic Klein-Gordon equation in three spatial dimensions. 



APPENDIX: PROOF OF EQ. ([76]) 

In this appendix we will prove Eq. (fZq): 

V^5o _ rfln(pj) 
p dl 

where J is the Jacobian of the transformation from the curvilinear coordinates 

I, a, P (which are used to describe the wave fronts S'o(r) = constant) to Cartesian 

coordinates. In order prove this relation, we must first prove the following 

lemma: 

dJ ,„ /VS'o\ ,„ dx 



which states that the logarithmic derivative of the Jacobian J along a congru- 
ence of trajectories is equal to the divergence of the tangent vector field of the 
congruence. 

A standard result, proven in many relativity texts pQ], is that for any matrix 
Qij, with determinant a and inverse a*-' — A^''/a, where A^^ is the signed cofactor 
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of aij , we have 



da ■■ dai. 



dx'^ 



dx'^ 



,^3 ^«y 



Let 



J 



9x 



for aij symmetric. 



(100) 



(101) 



d{l,a,l3) 
be the determinant of the transformation matrix from the curvihnear ray co- 



ordinates (I, a, (3) to the Cartesian coordinates x. Writing out Eq. (100) with 



Jij = dx^ /dx'^ , with a;* as Cartesian coordinates and a;" as curvihnear 



coordinates yields 



dJ dx'^ d ( dx'^ 



dx'^ 



dx"^ dx'^ \dx'^ 



= J 



dx'^ d 



dx^ dx'i \ dx''^ 



dx^ 



-'i?{^^- c"^' 



where in the second step we have switched the order of partial differentiation, 
and in the third step we have used the chain rule. If we set k — \ with x^ = I, 
the expression in the last parentheses above is dx/dl, in component form. Thus, 
with k = 1, Eq. (102) is just Eq. (pS) in component form. 

Now, using Eq. (p9|), we want to prove Eq. ([Tq). We begin with the first 
equality in Eq. (|99|) and expand out the divergence to obtain 



din J 
dl 



VS'o Vp 
P P 



V^So 



Using the ray equation, Eq. (69), in the form 

dx VS'o(r) 



we get 



dl 

din J dx 

dl " ^~dl 



Vlnp 



V^^o 



or, with use of Eq. (|70|), 



din J 
dl 



dhip 
dl 



(103) 



(104) 



(105) 



(106) 



Solving for V Sq/p in Eq. (106) yields the desired result, Eq. (|76|). 
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